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ABSTRACT 


This  thesis  suggests  a  method  to  estimate  the  current 
value  of  an  ionospheric  parameter.  The  proposed  method  is 
based  on  the  known  variability  of  the  observed  current  values 
near  path  and  utilizes  data  derived  from  ionospheric  sampling 
measurements.  Analysis  of  errors  is  provided  in  Single-Site- 
Location  High-Frequency  Direction  Finding  (SSL-HFDF) ,  arising 
from  ionospheric  irregularities  such  as  Es  (sporadic  E)  , 
ionospheric  tilts,  and  traveling  ionospheric  disturbances 
(TIDs)  .  The  characteristics  of  Es,  tilts  and  TIDs  for  mid¬ 
latitudes  are  summarized  in  tables.  The  spatial  and  temporal 
coherence  of  ionospheric  variedoilities  and  irregularities  is 
analyzed  over  the  electron  density.  Practical  results, 
measurements,  and  studies  are  presented  on  SSL-HFDF.  A  survey 
of  characteristics  of  the  ionosphere  in  the  equatorial  region 
is  also  provided.  Finally,  some  recommendations  are  given  to 
maximize  the  applicability  of  the  proposed  method. 
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I.  INTRODUCTION 


A.  MULTIPLE  SITE  HIGH  FREQUENCY  DIRECTION  FINDING 

Interception  of  high  frequency  (HF)  conununications  is  one 
means  of  gaining  additional  intelligence  in  a  military 
engagement.  The  interception  of  HF  communications  can  also  be 
used  to  locate  the  emitter,  the  knowledge  of  which  may  be  of 
strategic  or  tactical  value.  Early  radio  direction  finding 
employed  geometrical  triangulation.  Geometrical  triangulation 
involves  the  measurement  of  the  bearing  or  azimuth  angle  cf 
arrival  (measured  clockwise  from  the  geographic  north)  of  the 
target  signal. 

A  network  of  at  least  two  bearing  measuring  stations  is 
often  used  to  locate  HF  emitters.  Position  is  usually  derived 
from  the  intersection  of  two  or  more  bearing  estimates  which 
are  assumed  to  be  straight  line  ray  paths.  Fixing  accuracy  is 
increased  by  using  more  than  two  direction  finding  (DF) 
stations  and  indicates  how  closely  the  intersections  of  the 
various  rays  coincide.  The  most  likely  source  location  is 
derived  by  a  geometrical  center-of-gravity  (see  Figure  1) .  In 
Figure  1,  "FIX  n/m"  means  fix  using  the  DF  stations  DFn  and 
DFm. 
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Figure  1  Multiple  Site  Location  concept 

The  validity  of  the  straight-line  path  assxomption  as  well 
as  the  accuracy  of  the  azimuthal  angle-of-arrival  (AOA) 
measurements  regulate  the  degree  of  coincidence  of  various  ray 
intersections  and  hence  the  fixing  accuracy. 

At  frequencies  between  3  and  30  MHz,  over  paths  extending 
well  beyond  the  optical  horizon,  the  straight-line  ray  path 
assumption  is  not  always  valid.  The  most  likely  mode  of 
propagation  in  this  situation  is  typically  the  skywave  or 
ionospherically  propagated  wave. 

When  assuming  horizontal  stratification  of  the  refracting 
portion  of  the  ionosphere,  the  ray  will  bend  back  toward  the 
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earth  along  a  path  within  the  plane  of  incidence  and  the 
straight-line  ray  path  remains  valid.  Unfortunately  horizontal 
stratification  is  only  a  first-order  approximation  of  the 
ionosphere.  Large  ionization  gradients  are  present  because  of 
the  direct  effect  of  solar  radiation  on  the  ionosphere,  and/or 
ionospheric  disturbances  of  various  scales  and  magnitudes. 
Lateral  deviations  in  the  ray  path  occur  because  of  ionization 
gradients  which  have  components  transverse  to  the  ray  path, 
invalidating  the  straight-line  ray  path  assiimption.  It  seems 
evident  that  substantial  fixing  errors  remaining  in  the  HF 
geometrical  triangulation  procedure  are  due  largely  to  lateral 
deviations  in  the  ionospherically  propagated  ray  path. 

B.  SINGLE  SITE  HIGH  FHEQUENCY  DIBECTION  FINDING 

The  predominant  method  for  radio  wave  propagation  at  HF  is 
ionospheric  reflection.  The  use  of  SSL-HFDF  techniques 
requires  compensation  for  the  lateral  deviation  error  inherent 
in  the  triangulation  of  skywave  signals  and  to  avoid  the 
tactical  problem  of  installing  and  coordinating  angle-of- 
arrival  measurements  at  two  or  more  widely  separated  DF 
stations . 

SSL-HFDF  is  one  technique  employed  to  solve  the  problem  of 
determinating  the  location  (fixing)  of  HF  emitters.  As  the 
SSL-HFDF  station  makes  angle  measurements  in  the  vertical  and 
horizontal  plane  (see  Figure  2),  the  location  of  the  emitter 
is  evaluated  by  the  combination  of:  (1)  The  measured  bearing 


3 


<azimuthal  angle  measurement  that  was  previously  required  in 
the  triangulation  techniques);  (2)  The  observed  angle-of- 
elevation  (the  angle  made  by  the  arriving  ray  with  the  ground 
at  the  SSL-HFDF  station) ;  (3)  Ionospheric  models  available; 
and  (4)  Raytracing  techniques. 


Figure  2  Single  Site  Location  concept 

The  emitter's  location  is  obtained  by  measuring  the 

azimuth  bearing  and  the  range  is  calculated  by  tracing  a  ray 

«  / 

at  the  observed  angle  of  elevation  bade  towards  the  emitter, 
using  appropriate  ionospheric  models  and  raytracing  techniques 
[Ref.  1]. 


Any  emitting  HF  source  yields  two  types  of  waves:  one 
which  propagates  along  the  earth-air  interface  (the  surface 
wave  or  groundwave)  and  the  other  which  propagates  through  the 


atmosphere  (the  direct  wave) ,  eventually  to  be  reflected  off 
the  ionosphere  (the  skywave) .  The  proportion  of  the  one  wave 
to  the  other  for  any  given  emitter  depends  on  antenna  design, 
orientation,  distance  and  ground  conditions.  The  groundwave  is 
the  strongest  signal  received,  when  the  distance  is  20  km  to 
100  km,  depending  on  the  topography,  soil  type,  and  soil 
conductivity.  In  this  region  the  reflected  skywave  is  of  near¬ 
vertical  incidence  and  yields  little  useful  range  or  direction 
information.  For  short-distance  HFDF,  a  system  utilizing  the 
groundwave  would  be  the  most  applicable.  At  ranges  greater 
than  100  km,  and  most  optimally  at  ranges  of  a  few  hundred 
kilometers,  the  ground-based  HFDF  system  using  skywave  should 
be  employed.  Skywave  propagation  is  fundamental  to 
understanding  the  SSL-HFDF  process  [Ref.  2j . 

C.  THE  ROLE  OF  THE  IONOSPHERE 

The  ability  to  trace  the  ray  through  the  ionosphere 
depends  on  the  accuracy  of  the  model  of  the  refracting  portion 
of  the  ionosphere.  In  particular,  knowledge  of  the  vertical 
electron  density  profile  in  the  proximity  of  the  refracting 
ray  is  imperative. 

Distribution  of  electron  density  in  the  ionosphere  is 
strongly  affected  by  the  presence  of  horizontal  gradients  that 
cause  more  problems  than  in  the  case  of  the  triangulation. 
Accuracy  of  a  resulting  fix  by  the  SSL-HFDF  system  is  directly 
affected  by  the  longitudinal  as  well  as  the  transverse 


5 


components  of  horizontal  ionization  gradients .  Transverse 
components  affect  the  estimate  of  the  bearing  of  the  target 
from  the  SSL-HFDF  site,  while  the  longitudinal  components 
affect  the  estimate  of  the  range  of  the  target  from  the  site. 
The  backward  ray-tracing  procedure  must  rely  upon  incomplete 
information  regarding  the  current  status  of  the  ionosphere  as 
a  consequence  of  spatial  and  temporal  ionospheric  variations. 
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II.  REFRACTION  POINT  ESTIMATION  ANALYSIS 


A.  GENERAL 

For  ideal  performance  of  the  SSL-HFDF  station,  near-real- 
time  information  is  required  on  the  state  of  the  ionosphere. 
Ionospheric  sounding  is  most  frequently  used  to  accomplish 
this  in  the  region  near  the  path.  Quantifying  the  spatial  and 
temporal  irregularities  of  the  ionosphere  becomes  very 
important.  Equally  important  is  the  estimation  of  spatial  and 
temporal  ranges,  which  allow  ionospheric  sounding  information 
gathered  at  one  point  to  be  extrapolated  to  another  point  with 
minimal  loss  of  position  location  accuracy. 

The  fundeunental  question  is,  "If  the  state  of  the 
ionosphere  can  be  determined  at  one  particular  point,  over 
what  spatial  ranges  can  that  information  be  transferred,  and 
for  what  time  period  is  it  valid  ?"  [Ref.  €] . 

The  principal  aim  of  this  chapter  is  to  present,  in 
theory,  one  process  for  the  estimation  of  the  current  value  of 
an  ionospheric  pareuneter.  The  applicad>ility  of  this  process 
should  be  investigated  by  using  real  ionospheric  data 
measurements,  which  were  not  avail2d>le  for  this  research. 

The  limitations  concerning  spatial  ranges  and  time  period 
will  be  discussed  in  the  final  chapter. 
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B.  THE  BSTIM21TI0N  PROCESS 


1 .  Introduction 

The  problem  considered  here  is  that  of  estimating  the 
current  value  of  an  ionospheric  parameter  on  the  basis  of 
observed  current  values  of  the  parameter  at  other  locations . 

In  this  study  the  following  notations  and  assumptions 
will  be  made: 


•  The  location  for  which  the  estimation  is  to  be  performed 
is  given  the  label  0; 

•  The  locations  for  which  current  observed  values  are 
avail2±)le  are  designated  l,2,3,...n; 

•  The  expected  value  (mean)  and  the  variance  of  the 
parameter  at  each  location  0,1, 2, 3, . . .,n  are  assumed  to  be 
known; 

•  The  covariance  of  the  parauneter  between  each  pair  of 
locations  is  assumed  to  be  known. 


The  way  to  estimate  the  previous  parauneters  can  be 
summarized  as  follows: 


•  The  mean  value  at  location  0  may  be  estimated  in  term  of 
the  median  using  the  availadDle  ionospheric  predictions 
published  by  specialized  organizations  around  the  world; 

•  The  variance  at  location  0  is  estimated  from  the 
geographical  dependence  observed  in  the  variance  at  the 
other  locations; 

•  The  covariances  involving  0  are  estimated  from  the 
distance  dependencies  [Ref.  3] . 
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2 .  Expected  Values  and  Variances 

Estimation  of  the  deviation  of  a  particular  parameter 
from  its  mean  value  at  location  0  is  as  described  in  the 
following  analysis  [Ref.  3] . 

It  is  necessary  to  calculate  the  predicted  current 
deviation  Xj,  at  0,  which  is  assumed  to  be  a  linear  function  of 
the  observed  current  deviations  at  the  n  other  locations, 
i  .  e . , 

n 

Xp  =  n'U  g 

where 

are  the  observed  deviations, 
ki  are  the  prediction  coefficients,  and 
n  is  the  number  of  the  location. 

The  goal  is  to  minimize  D,  the  expected  mean  square 
prediction  error  at  0,  which  is  given  by 

D=  (Xp'Xo)2  ,  (2) 

where  Xq  represents  the  true  deviation  at  0. 
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Then 


D  =  [  (  )  -Xq  ]  2 


(3) 


Where  the  bar  over  the  expression  represents  the  mean 
or  expected  value. 

To  minimize  D,  the  partial  derivatives  with  respect  to 
the  n  ki's  must  be  zero.  Thus,  the  solution  of  the  following 
system  of  n  equations  yields  the  ideal  values  for  the 

-^=2  [(  ^  ~^0  1  (  n'^X^  )  =  0  ,  (^) 


-^=2  [(  n'^'^k^x^  ) -Xq  ]  (  n‘^X2  )  =  0  ,  (5) 


and 


[  (  k^x^  )  -Xj,  ]  (  n''^x„  )  =  0 


(6) 


m.  C*am  1  (Onm  location) 

When  estimating  a  parameter  at  location  0  with 
only  one  location  (n=l) ; 
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(7) 


-  XqX^  =  0 


or 


XqXj  _ 


XiXi 


(8) 


If 


S^J  ~  X^Xj 


then 


= 


(9) 


(10) 


where 

S(,i  is  the  covariance  between  locations  0  and  1,  and 
Sii  is  the  variance  at  location  1. 


The  constant  can  also  be  described  as  a  function 
of  the  correlation  coefficient  r^i. 


11 


So  from 


r 


01 


Cov  (0,1) 


(11) 


Cov  (0,1) 


^01 


•^01®0®l 


(12) 


and  hence 


Sii  = 


(13) 


The  coefficient  can  be  defined  as 


(14) 


where  the  correlation  coefficient  rgi  between  locations  0  and 
1  [Ref.  4],  is  given  by 


r 


01 


COV^(0. 1) 


(15) 


If  Soof  the  variance  at  0,  is  equal  to  s^,  then 
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■^01 


®01  _  ®01 


K  =  .  (17) 

Hence  the  mean  square  error  D  can  be  calculated  as 

D  =  [  k^x^-Xo  =  [  r^^x^-Xo  =  roiX^'2-i'oi-^i'^o^V 


where : 


so  that 


2  -^1  X,Xr, 

-  2ro,-i^  -  1 
Xq  Xo 


_ 


*^1^0  .  Sio  _  Soi  _ 


D  =  Xo  [  roi  -  2roi  +  1  ]  =  Xo  [  1  -  Toil 


Sqo  (1  ^01^^ 
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Finally, 


where 

b. 

dP 

dk^ 

and 

dP  ^ 
dk^ 

We  have 

and 


P  =  o,2(l-r„,2) 


(24) 


is  the  variance  at  location  1 . 


Camm  2  (Two  locations) 

When  estimating  0  with  two  locations  (n=2) : 


2  I  2  -1  _  _  _ 

-  =  [k^x^^+k2X^X2)  -^0-^1 

2  2 


(25) 


2  [  —  ( k-^x.^  ^iCjXj )  ~Xq]  X2 


=  [k^x^X2+k2X2^)  -x^X2 


(26) 


dP 

dk^ 


1/2  -Sqi  =  0 


(27) 


dP 

dk2 


\  /  “Z  ^  k  ■^S-^2^  ^2^22^  ®02  ® 


(28) 
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so  that 


= 


■  12 


'11  “^12 


(1--^12^) 


(29) 


K  = 


®02  X  _  t  ^01  X  ,  2 

12 


(zsi)-(s:2i)r 

^22  ^12 


(30) 


D  = 


^  -Xq 


(31) 


hence, 


—  ^ ^k^k^x^x^^x^^ -k-^x^x^-k^x^x^  ,  (32) 


^  4  (^1^Sii'*’^2^®22)  ■*■  *2  ^1^2^12  ■*■^00  “^1^01  ~^2®02  •  (33) 

Solution  of  this  system  when  n=2  can  be  simplified 
by  assuming  eq[ual  variances. 


®11  ®22  ~  ®00  *  ®XX 


(34) 


and  equal  covariances  with  respect  to  0, 


®01  ~  ®02  ~  ®0x 


(35) 
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so  that 


=  ^2  =  = 


2  s, 


Ox 


(36) 


and 


D  =  ^s^-2k^So^ 

which,  substituting  for  k,,  gives 


(37) 


D  =  2- 


^QX 


+  S  -4  * 

«rw  ^ 


^Ox 


'  (s,,+s,,) 


"■xx  ■^12' 


(38) 


(it-Ja) 

®xx 


1  —  r  2 
-^Ox 


2 

l+ri2) 

(39) 


This  formula  reduces  to  the  case  n=l  when  ri2=l . 

c.  Cmmm  3  (Multxplm  locations) 

When  the  prediction  is  based  upon  observations  at 
more  than  one  location,  the  general  approach  to  the  expected 
mean  square  prediction  error  (See  Appendix  A)  is: 
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D  =  0,2  (l-i?o„2) 


(40) 


where  is  the  coefficient  of  multiple  correlation  of  x, 

with  relation  to  {  x^,  Xj,  Xj  }. 

C.  SUMMARY  OF  THE  ESTIMATION  PROCESS 


1st  STEP;  Calculate  the  variance  (Oi^)  of  the  par8uneter  at 
each  location  i=l,2,...,n  given  by  the  formula 


_  i-i 

Oi  -  - 


_  _ 


Lx. 

iill— 
n 


(41) 


n-1 


n-1 


where  n  *  number  of  measurements  at  location  i, 
x^  s:  measurement  at  location  i,  and 
Hi  =  mean  value  of  the  observations  at  location  i 

[Ref.  4]. 

2nd  STEP;  Calculate  the  variance  {<5^)  of  the  parauneter  at 
location  0  as  in  the  first  step  or,  when  no  measurements  are 
availa±)le,  estimate  its  value  from  the  spatial  (geographical) 
interpolation  of  the  variances  at  locations  i=l,2,...,n. 
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3rd  STBP ;  Calculate  the  covariances  and  the  correlation 


coefficients  for  each  pair  of  locations  including  location  0. 

The  correlation  coefficient  (r^y)  is  the  simple  linear 
correlation  coefficient  given  by  the  formula  [Ref.  5] 


= 


(  "  1 

(  n  \( 

17 

- 

[ill _ 

A 

n 


n 

Vi-I 


►  X 

V-rf  } 


(42) 


Vi-i 


/  n 

U'l  - 


The  covariance  (COV(X,Y))  for  each  pair  of  locations  is 
given  by  the  formula  [Ref.  4] 


COV{K,  Y)  = 


(43) 


4th  3TKP :  Solve  the  correlation  matrix  (see  Appendix  A) 
and  calculate  the  prediction  coefficients  (k^)  and  the 
coefficient  of  multiple  correlation  (Ro*^) 


kj  = 


(j  =  1,2 . n) 


(44) 


R 


2 

On 


1- 


_l^ 

^00^00 


(45) 
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5th  STEP :  Calculate  the  current  deviation  (x^)  of  the 


parameter  for  each  location  i=l,2,...,n. 

where  =  deviation  of  the  ith  meacurement  from  the  mean 
(ith  deviation); 

m^  =  measured  or  observed  value  at  location  i; 

|ii  =  actual  mean  value  at  location  i. 


6th  STEP ;  Calculate  the  estimated  deviation  (x^)  of  the 
parameter  at  location  0: 


j[*i 


(47) 


7th  STEP :  Calculate  the  predicted  value  (P)  of  the 
parameter : 


P  = 


(48) 
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Bth  STKP ;  Calculate  the  expected  mean  square  prediction 
error  (D)  for  one  station  and  for  multiple  stations. 

For  one  location,  the  error  is 

£)  =  Oo(l-roi)  , 

and  for  multiple  locations  it  is: 

D  =  ,  (50) 

where  r<,i  =  Correlation  coefficient  between  locations  0 

and  1 ,  and 

Ron  =  Multiple  correlation  coefficient  among 
location  0  and  locations  l,2,...,n. 
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III.  ANAI.YSIS  OF  ERRORS  IN  SSL-HFDF 


A.  OVERVIEW 

Current  SSL-HFDF  systems  occasionally  produce  errors  of 
tens  of  kilometers  or  worse.  These  errors  arise  from  three 
main  areas:  ionospheric  variability  and  irregularity,  locator 
system  size  limitations,  and  problems  with  data  acquisition, 
processing  and  interpretation.  Of  these  areas  the  ionosphere 
is  the  single  largest  source  of  error. 

Measurement  of  the  azimuthal  and  elevation  angles  of  the 
incoming  signal  and  determination  of  the  height  of  the 
ionospheric  reflecting  layer  are,  in  principle,  the  procedures 
used  to  locate  the  HF  emitter. 

In  practice,  there  are  several  types  of  ionospheric 
irregularities  which  distort  the  otherwise  straightforward 
picture  of  a  uniform,  concentric,  smoothly-reflecting 
ionosphere.  Three  phenomena  can  be  pointed  out  because  of 
their  significance  among  various  causes  of  irregularities  and 
their  effects  on  SSL-HFDF;  sporadic  E  (Es) ,  ionospheric  tilts, 
and  traveling  ionospheric  disturbances  (TIDs) . 

1 .  Sporadic  E  (Es) 

Because  of  the  thickness  of  the  layers  and  the 
sharpness  of  the  vertical  gradients.  Sporadic  E  is  often  an 
aid  in  communication  and  often  provides  a  better  reflection 
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surface  than  the  F  region.  Unfortunately,  the  inconvenience  of 
Es  is  that  it  introduces  uncertainty  as  to  which  ionospheric 
layer,  E  or  F,  is  reflecting  the  signal  and  provides 
additional  opportunities  for  multimode  reflections  of  signals. 
Signals  from  the  source  may  suffer  reflections  from  both  E  and 
F  layers,  thus  causing  additional  problems  with  resolution  or 
modes . 

The  properties  of  Es  are  summarized  in  Table  I 
[Ref.  6]  for  the  mid-latitudes . 

TABLE  I  -  PROPERTIES  OF  Es  FOR  THE  MID -LATITUDES 


Structure 

(1)  Patches  of  enhanced  electron  and  ion 
density  often  hundreds  of  km  in  horizontal 
extent . 

(2)  Vertical  thickness  generally  1-2  km  at 
an  altitude  of  100-110  km. 

(3)  Patches  of  high  density  often  embedded 
in  larger,  lower  density  patches. 

Origin 

Wind  shears,  probeQ:>ly  from  propagating 
acoustic-gravity  waves. 

Motion 

Patches  (not  plasma)  generally  move  50- 
100  m/s,  no  preferred  direction. 

Duration 

Several  minutes  to  several  hours . 

Occurrence 

(1)  More  frequent  during  day,  with  peak 
occurrence  before  noon,  and  in  some 
locations  a  secondary  peak  near  sunsets 
likely  in  summer. 

(2)  Frequency  of  occurrence  (detection) 
more  likely  for  lower  radio  frequencies . 

2 .  Ionospheric  Tilts 

Ionospheric  tilts  refer  to  any  deviation  from  the 
horizontal  plane  of  the  contours  of  constant  electron  density. 
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whether  caused  by  large-scale  phenomena  such  as  solar 
ionization,  or  due  to  more  transient  and  localized  TIOs. 

Miscalculation  of  the  virtual  height  of  the  reflection 
point,  and  error  in  estimating  the  angle  of  arrival,  in 
particular  the  elevation  angle,  are  the  dominant  effects  of 
ionospheric  tilts. 


Figure  3  Effect  of  a  Longitudinal  Tilt  component 

The  diurnally-varying  solar  ionization  rate  in  the  F 
region  causes  the  most  regular  and  predictedjle  ionospheric 
tilt.  This  is  most  observzible  at  sunrise  and  sunset,  so  it  can 
be  predicted  and  described  [Ref.  6] . 
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The  properties  of  ionospheric  tilts  are  suromarized  in 


Table  II  [Ref.  6]  for  the  mid-latitudes. 


TABLE  II  -  PROPERTIES  OF  IONOSPHERIC  TILTS  FOR  THE 
MID -LATITUDES 


Structure 

(1)  Horizontal  gradient  in  electron 
density  over  distances  on  the  order  of 

1000  km. 

(2)  Found  at  F  region  altitudes  (>140  km) . 

Origin 

Daily  variation  of  solar  ionizing 
radiation . 

Motion 

Pattern  shifts  with  sun's  diurnal  motion. 

Duration 

On  the  order  of  1-2  hours. 

Occurrence 

Daily  near  sunrise  and  sunset. 

3.  Traveling  Ionospheric  Disturbances  (TIDs) 

Errors  in  the  measured  azimuthal  and  elevation  angles 
(angles  of  arrival) ,  and  in  the  virtual  height  of  reflection 
are  the  main  effects  of  TIDs. 

The  spectrum  of  TIDs  can  be  placed  in  at  least  three 
distinct  categories;  large-scale,  medium-scale  and  small-scale 
[Ref.  6] . 


•  Large-scale  TIDs,  according  to  acoustic  wave  theory,  are 
associated  with  a  discrete  spectrum  of  guided  waves  whose 
modes  are  excited  only  by  upper  atmospheric  sources  and 
whose  horizontal  speeds  are  s\ibstantially  greater  than  the 
(lower  atmospheric)  speed  of  sound. 

•  Medium-scale  TIDs  are  associated  with  a  spectriim  of 
freely-propagating  internal  waves  which  can  be  excited  by 
sources  at  any  altitude  and  whose  horizontal  speeds  are 
less  than  the  speed  of  the  sound.  Medium-scale  TIDs  are 
much  nore  common . 
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•  Small-scale  TIDs  are  more  likely  the  extension  to  higher 
frequencies  and  smaller  size  of  the  medium-scale  TIDs. 
They  are  generally  below  the  Fresnel-zone  size  of 
ionospheric  sounders  and  thus  have  not  been  as  well 
documented. 


Figure  4 


Traveling  Ionospheric  Disturbance  (TID) 


TaJDle  III  [Ref.  6]  svunmarizes  the  properties  of  these 
different  categories,  and  Ted>le  IV  [Ref.  6]  gives  the 
magnitude  of  the  errors  in  position  location  which  the 


ionosphere  can  cause  for  selected  ranges .  A  quick 
thumb"  seems  to  be  10  km  or  10%  of  range,  whichever 
[Ref.  6]  . 


"rule  of 
is  worse 
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TABLE  III  -  PROPERTIES  OF  TIDs 


TID 

LARGE-SCALE 

MEDIUM-SCALE 

SMALL-SCALE 

Wavelength 

(1)  >1000kin 

(1)  lOs  to 

(1)  <10  km 

and 

horizontal 

lOOs  km 

horizontal 

structure 

wavelength 

horizontal 

wavelength 

wavelength 

(2)  Wave- 

(2)  Wave- 

(2) Structure 

front  width 

front  width 

not  well 

near  1000 

100s  to  1000 

resolved 

km 

km 

(3)  Phase 

(3)  Phase 

fronts 

fronts 

tilted 

tilted  30“- 

nearly 

60®  from 

horizontal 

vertical 

(4)  Retains 

(4)  Do  not 

shapes  over 

retain 

1000s  km 

shapes  at 
d>100  km 

Motion 

>  300  m/ s 

100-250  m/s 

100-250  m/s 

North  to 

variable 

(estimated) 

South 

directions 

variable 

with 

seasonal 

trends 

directions 

Period 

30  min. -3 

10-100  min. , 

<10  min. 

hr. 

several 

long,  trains 

usually  1-3 

cycles  to 

to  families 

cycles 

trains 

as 

wavelength 

decreases 

Frequency 

Infrequent, 

Daily,  more 

Daily 

of 

less  than 

common 

Occurrence 

daily 

daytime 

Source 

(1)  Events 

(1)  Tropos- 

(1)  Probably 

in  the 

pheric 

tropospheric 

auroral 

phenomena 

(not  well 

zone 

established) 

(2)  Strong 

(2)  Upper 

correlation 

atmospheric 

with 

&  polar 

magnetic 

winter 

activity 

sources 
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TABLE  IV  -  MAGNITUDE  OF  ERRORS  IN  HFDF  DUE  TO  ERRORS  IN  AOA 


FOR  E  AND  F  LAYERS 


ACTUAL  RANGE  (km) 

200 

300 

HEIGHT  OF  REFLECTING  LAYER  (km) 

105 

250 

105 

250 

1  Range  Error  (km)  for:  | 

1°  elevation  angle  uncertainty 

6.9 

9.8 

10.9 

11.5 

3°  elevation  angle  uncertainty 

20.6 

29.3 

32.7 

34.4 

1  Cross-Range  (Azimuthal)  Error  (km)  for:  | 

1°  azimuthal  angle  uncertainty 

3.5 

3.5 

5.2 

5.2 

3°  azimuthal  angle  uncertainty 

10.5 

10.5 

15.6 

15.6 

1  Range  Error  (km)  for:  | 

1  10  km  height  uncertainty 

27.6 

11.5 

B.  SPATIAL  AMD  TEMPORAL  COHERENCE  OF  IONOSPHERIC 
VARIABILITIES  AND  IRREGULARITIES 

1 .  Electron  Density :  Tcnporal  and  Spatial  Variabilitiaa 
The  ionosphere  displays  variation  in  electron  density 
on  many  temporal  scales.  The  11-year  solar  cycle  dependency  in 
addition  with  the  seasonal  and  monthly  changes  in  electron 
density  are  the  most  useful  scale  references. 

The  day-by-day  variability  is  the  most  difficult  to 
model.  As  an  example,  some  studies  show  that  the  variad^ility 
of  daily  hourly  values  of  foE  and  foFl  (see  Glossary)  about 
the  monthly  hourly  values  expressed  in  terms  of  the  standard 
deviation  is  on  the  order  of  5  to  10  percent.  The  saune  study 
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provides  evidence  that  the  daily  variability  of  foF2  (see 
Glossary)  is  between  10  and  20  percent  [Ref.  7] . 

The  equatorial  and  polar  regions  of  the  globe  show  a 
daily  variaUbility  stronger  than  that  in  mid-latitude 
locations.  This  phenomena  is  attributed  to  the  daily 
variability  observed  in  processes  such  as  electrodynamic 
drift,  diffusion,  neutral-air  winds,  and  particle 
precipitation  that  compete  for  control  of  the  ionization 
distribution  in  the  low-  and  high-latitude  F  region. 

It  is  very  important  to  emphasize  that  the  height  of 
the  maximum  electron  density  (hm)  in  the  different  ionospheric 
regions  also  varies  with  time  and  location,  and  the  hmF2  tends 
to  show  the  most  significant  changes  of  all  the  ionospheric 
maximum  heights.  At  mid- latitudes,  hmF2  is  usually  on  the 
order  of  300  to  400  km  during  nighttime  hours  and  between  250 
and  350  km  during  daytime  hours,  whereas  at  equatorial 
latitudes,  the  daytime  values  are  100  to  150  km  higher  than 
the  nighttime  values. 

Regardly  spatial  variability,  the  electron  density  is 
dependent  on  the  solar  zenith  angle,  as  a  general  rule,  and 
its  maximum  latitudinal  distribution  is  a  maximum  near  the 
sub-polar  point.  The  F  region  amd  E  region  (Es)  have  a 
significant  latitudinal  dependence,  and  in  the  F2  region,  the 
electron  density  displays  complicated  latitudinal  behavior 
[Ref.  8] . 
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2 .  Blttctron  Density  Irregularities 

The  analysis  of  the  irregular  electron  density 
structures  in  the  ionosphere  in  the  F  and  E  regions  are 
summarized  in  TaUole  V  [Ref.  8]  . 


TABLE  V  -  IONOSPHERIC  ELECTRON  DENSITY  IRREGULARITIES 


REGION 

LATITUDE 

IRREGULARITIES  IN  THE  ELECTRON  DENSITY 

LOW 

.  In  evening  hours,  continuing 
throughout  night  hours 
.  Tends  to  decrease  during  times  of 
geomagnetic  disturbances 
.  Significant  longitudinal  variadsility 
as  well  latitudinal  and  temporal 
dependencies 

.  Causes  fading  on  HF  circuits 

F 

MIDDLE 

.  Associated  with  Spread  F  phenomenon 
.  Low  magnitude,  duration,  and  effect 
on  communications  system  performance 

HIGH 

.  It  appears  more  or  less  routinely 
.  More  intense  during  hours  of 
darkness 

.  Greater  than  at  low  latitudes 
.  It  causes  fading  on  HF  circuits 

LOW 

.  Generally  are  characterized  as  Es 
.  Present  during  high  %  of  daytime 
.  Latitude  of  maximum  occurrence 
varies  with  longitude 

E 

MIDDLE 

.  Are  generally  characterized  as  Es, 
most  frequently  in  summer  daytime 
.  The  diurnal  variation  has  maxima  in 
the  mid-morning  hours  and  near  sunset 

HIGH 

.  Irregularities  are  quite  common 
.  Associated  with  auroral  activity 
are  mainly  a  nighttime  phenomenon,  and 
when  occurring  in  thick  layers  they 
produces  radio  wave  retardation 
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3.  Practical  Results  on  SSL-HFDF 


The  reflection  point  of  the  intercepted  HF  signal  is 
far  from  the  SSL-HFDF  station.  For  a  spatially  and  temporally 
uniform,  or  at  least  slowly  varying  ionosphere,  this  would 
introduce  no  problem.  In  reality,  however,  ionospheric 
irregularities  reduce  the  usefulness  of  information  gathered 
at  one  point  when  transferred  to  another  point. 

By  cross-correlating  the  angle-of-arrival  de"'  ^ions 
of  the  signals  received  from  pairs  of  geographically  distinct 
emitters,  and  assuming  a  one-hop  propagation  path,  some 
practical  results  can  be  obtained.  When  the  variations  in  the 
angle  of  arrival  at  one  location  are  reproduced  at  the  second 
location  "t"  seconds  later,  the  cross-correlation  function  is 
maximized.  Interference  between  signals  from  different 
sources,  or  from  the  szune  source  travelling  different  paths, 
or  due  to  a  single  frequency  component  of  the  composite 
disturbance  can  result  in  the  observed  decorrelation. 

Good  spatial  and  temporal  coherence  are  observed  from 
ionospheric  tilts  due  to  solar  influence  or  from  large-scale 
TIDs,  maintaining  their  shapes  over  long  distances  and  for 
times  on  the  order  of  one  hour  or  more.  The  spectrum  of 
medium-scale  TIDs  is  superimposed  on  these  migrating  regular 
waves .  The  most  frequently  occurring  ionospheric 
irregularities  and  the  ones  showing  the  least  spatial  and 
temporal  coherence  are  the  medium-scale  TIDs,  and  these  seem 
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to  place  the  most  rigorous  limitations  on  accuracy  when 
extrapolating  the  state  of  the  ionosphere. 

Sporadic  E  can  be  analyzed  under  the  concepts  of  spatial 
and  temporal  coherence,  or  decorrelation,  only  when  the  patch 
length  and  the  flux  of  the  Sporadic  E  layer  will  give  some 
estimate  of  how  long  the  phenomenon  is  expected  to  remain  at 
any  one  given  point . 

Table  VI  [Ref.  6]  provides  a  summary  of  this  group  of 
properties  for  ionospheric  irregularities. 
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TABLE  VI 


SUMMARY  OF  IONOSPHERIC  IRREGULARITIES 


TYPE 

SPATIAL 

DECORRELATION 

TEMPORAL 

DECORRELATION 

MAJOR 

EFFECTS 

Es 

~l00s  km, 
depending  on 
the  size  of 
the  patch 

Minutes  to 
hours , depend . 
on  relative 
location  of 
patch  &  its 
drift  veloc. 

-  Uncertainty 
in  height  of 
reflecting 
layer 

-  Multimode 
propagation 

Ionospheric 

Tilt 

100s  to 

1000s  km 

Hour  or 
longer 

-  Uncertainty 
in  AOA  and 
height  of 
reflecting 
layer 

Large-scale 

TID 

1000s  km 

30  min  to 
several  hours 

-  Uncertainty 
in  AOA  and 
height  of 
reflecting 
layer 

Medium- 

scale 

TID 

50-100  km 

~5  minutes 

-  Uncertainty 
in  AOA  and 
height  of 
reflecting 
layer 

-  Multimode 
propagation 

C.  EQUATORIAL  REGION  CONSIDERATIONS 

The  considerations  discussed  here  do  not  give  a  thorougn 
study  of  all  problems  of  the  equatorial  ionospheric 
environment,  a  very  intricate  region,  but  they  do  reflect  the 
currant  knov/ledgement  of  this  medixim  and  point  out  some  of  its 
characteristics : 

(1)  Large  ionospheric  F-region  tilts  occur  along  great 
circle  paths  leaving  temperate  regions  and  passing  through 
equatorial  regions  in  early  evening  hours.  These  tilts,  both 
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positive  and  negative,  range  in  angle  from  about  1°  to  6° 
[Ref.  9]; 

(2)  Very  pronounced  modifications  of  the  modes  of 
propagation  of  HF  signals  along  paths  traversing  the 
equatorial  regions  in  evening  hours  are  produced  by  quite 
small  tilts  in  the  ionospheric  F2  region  [Ref.  9]; 

(3)  The  more  important  of  the  transequatorial  tilt-mode 
propagation  paths  occur  at  low  elevation  angles.  A  high  degree 
of  ray- focusing  commonly  occurs  for  elevation  angles  below  8°. 
These  focused  modes  lead  to  very  ste±>le  backscatter  echo 
configurations,  such  as  are  commonly  observed  for  anomalous 
transequatorial  echoes  [Ref.  9] ; 

(4)  There  is  a  permanent  daytime  belt  of  sporadic  E  of 
high  critical  frequency,  which  breaks  up  and  disappears  near 
sunset  [Ref.  10]; 

(5)  As  the  distance  between  two  points  increases,  there  is 
generally  a  decreased  correlation  between  the  pareuneters  at 
the  two  points,  tending  to  become  negative  at  the  distance  of 
the  stations  near  the  magnetic  equator  [Ref.  11]; 

(6)  The  spread  F  is  more  likely  occurs  between  about  2100 
and  0100  local  time,  earlier  when  sunspots  are  at  maximum.  Its 
occurrence  is  greater  during  local  summer  than  local  winter. 
Its  occurrence  probability  decreases  with  sunspot  number  [Ref. 
16]  . 
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rv.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  appliceibility  of  the  proposed  method  for  estimating 
ionospheric  parauneters  must  be  investigated  xising  data 
measurements .  The  predicted  values  calculated  should  be 
compared  to  real  measurements.  For  each  particular  environment 
(low,  middle  and  high  latitudes)  the  periodic-.ty  of  the 
measurements  and  some  adjustments  in  the  model  will  certainly 
be  necessary  to  find  the  optimum  predicted  parameter.  (For 
example,  in  the  equatorial  region  the  predicted  values  should 
be  calculated  at  sunrise,  noon,  sunset  and  nighttime) . 

B .  RECOMMENDATIONS 

The  following  recommendations  were  found  in  technical 
reports  and  stem  from  several  years  of  exhaustive  research  by 
specialized  organizations.  They  relate  to  maximizing  the 
applicaJoility  and  making  the  proposed  procedure  operative  for 
the  estimation  of  ionospheric  parameters. 

(1)  The  validity  of  the  ionospheric  data  taken  at  one 
point  is  lost  when  extrapolated  over  distances  of  50-100  km  or 
at  periodicity  greater  than  five  minutes  [Ref.  6] . 
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(2)  Data  from  a  single  ionosonde  should  never  be  used  to 
model  a  large  area  of  the  ionosphere.  The  utilization  of  an 
integrated  network  of  ionosondes  is  recommended  using  both 
vertical  and  oblique  sounding  [Ref.  6] . 

(3)  Analysis  of  the  relationship  between  R„^  (R„  = 
coefficient  of  multiple  correlation)  and  n  (number  of 
measurements  or  observations)  for  particular  configurations 
indicates  that  the  maximum  value  of  R„  (minimum  value  of  the 
expected  mean  square  prediction  error  D)  tends  to  be  reached 
at  low  values  of  n  [Ref.  11]  .  This  recommendation  suggests  the 
investigation  of  different  ionosonde  configurations  to  fit 
diverse  tactical  environments. 

(•  '  Installation  of  the  ionosonde  at  the  anticipated  mid¬ 
point  of  the  propagation  path  is  suggested  to  optimize  the 
measurements.  To  find  the  best  location  for  the  SSL-HFDF 
system  and  ionosonde  station,  a  preselection  of  range  and 
detection  is  recommended  over  which  the  fixing  will  occur 
[Ref.  6] ; 

(5)  Measurement  at  different  observation  times  is 
recommended  to  provide  data  for  parameter  estimation  as  a 
function  of  time  in  addition  to  simultaneous  time  geographical 
interpolation  [Ref.  11]; 
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(6)  Implementation  of  an  ionosonde  network  sometimes 
becomes  impractical,  as  for  exeunple,  in  an  ocean  region  or 
when  the  utilization  of  the  output  data  may  be  not  possible  at 
near-real-time  due  to  large  networks.  The  suggested  solution 
in  these  cases  is  to  use  techniques  to  obtain  data  from 
oblique  ionograuns  by  using  oblique  bistatic  sounding  or 
oblique  backscatter  sounding  techniques  [Ref.  12] . 
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APPENDIX  A  -  GENERAL  SOLUTION  OF  THE  CORRELATION  MATRIX 


Let  S  represent  the  correlation  matrix  [Ref.  3] ; 
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or 


i?  =  Soot  l-i2on^  1  =  o^[  l-i?o„2  ] 


(57) 


where 

S  is  the  correlation  matrix  with  rows  (i)  and  columns  (j) 
varying  from  0  to  n; 

Cij  is  the  cofactor  of  s^j  in  |sl;  that  is,  the  ijth 
cofactor  of  S  is  obtained  by  taking  the  determinant  of  the 
ijth  minor  and  multiplying  it  by  (-l)^"^^  [Ref.  13]; 

is  the  ijth  minor  of  the  n  x  n  matrix  S,  and  is  defined 
as  the  (n-1)  x  (n-1)  matrix  obtained  from  S  by  deleting  the 
ith  row  and  Jth  column  of  S  [Ref.  13]; 

|s|  is  the  determinant  of  S; 

is  either  the  covariance  between  locations  i  and  j 
(when  i  does  not  equal  j)  ,  or  the  variance  at  location  i,  a*i, 
(when  i=j)  ; 

Ro(i,2 a)  r  or  R^,  is  the  coefficient  of  multiple 

correlation  of  x^,  with  [Xj,  Xj,  ...,  x„};  and 

D  is  the  expected  mean  square  prediction  error. 
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APPENDIX  B  -  SKYNAVB  PROPAGATION  CONCEPTS 


A  complex  reflection  mechanism  of  ionospheric  layers 
controls  the  propagation  of  skywaves.  As  shown  in  Figure  5, 
the  propagating  wave  can  be  thought  of  as  a  ray  reflecting 
alternately  from  the  ionosphere  and  the  ground.  The  wave 
generally  follows  several  different  paths  at  the  saune  time, 
and  this  combination  of  paths  is  of  major  importance  in  the 
transmission  of  the  signal's  energy  between  two  points.  At  any 
given  time,  the  nature  of  this  multipath  propagation  is  also 
highly  dependent  on  ionospheric  conditions.  Axtnough  it  is 
called  "reflection",  the  phenomenon  which  causes  radio  waves 
to  change  direction  is  not  true  reflection,  but  a  case  of 
strong  refraction  in  a  downward  direction.  As  the  magnitude  of 
refraction  depends  upon  the  frequency,  the  effect  of  the 
refraction  decreases  with  increasing  frequency.  Hence,  there 
is  critical  frequency  for  a  given  ionospheric  profile  and 
angle  of  incidence,  beyond  which  reflection  back  to  earth 
cannot  occur.  According  to  observed  radio  properties,  the 
ionosphere  has  been  subdivided  into  three  regions  or  layers, 
D,  E,  and  F,  increasing  in  altitude  and  in  electron  density. 

The  D  region  spans  the  altitude  range  50  to  90  km,  with 
electron  density  varying  according  to  the  solar  zenith  angle. 
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Figure  5  Types  of  Oblique  Ray  Paths 
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In  the  case  of  the  HF  spectrum,  the  D  region  acts  principally 
as  an  attenuator  of  radio  waves. 

The  E  region  ranges  from  90  to  130  km  in  height,  and  the 
electron  density  encompasses  the  so-called  "normal"  and 
"sporadic  E"  (E,)  layers.  The  normal  E  layer  has  a  strong 
solar  zenith  angle  dependence,  and  the  E.  has  a  diurnal 
variation  with  location  on  the  globe.  The  normal  E  acts  as  a 
reflector  of  HF  waves,  particularly  during  the  daylight  hours, 
and  the  electron  density  associated  with  E,  acts  by 
reflecting  and/or  scattering  HF  waves. 

The  F  region  is  a  thick  region  that  extends  upward  from 
130  km.  It  is  subdivided  into  FI  and  F2  layers  because  its 
lower  part  displays  a  different  variation  than  the  upper  part. 
The  FI  layer  electron  density  reaches  its  maximtun  near  local 
noon  and  during  the  summer.  At  night  and  during  winter  there 
is  no  distinction  between  FI  and  F2  layers.  The  F2  layer 
displays  the  greatest  electron  density,  zmd  its  electron 
density  departs  significantly  from  a  simple  solar  zenith  angle 
dependence,  since  it  is  strongly  influenced  by  neutral-air 
winds,  electrodynamic  drifts,  and  diffusion  processes,  the 
maximiam  electron  density  tends  to  occur  after  local  noon  and 
its  structure  depends  upon  geomagnetic  latitude.  Long-distance 
HF  circuits  rely  primarily  upon  reflection  of  radio  waves  from 
the  F  region  during  both  day  and  night  hours. 
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The  altitude  and  density  of  the  F-region  peak  determines 
the  critical  frequency  for  propagation.  Especially  where 
sharp  gradients  occur  in  the  electron  density  profile, 
reflections  at  some  frequencies  can  occur  at  altitudes  well 
below  this  region  peak. 

The  interaction  of  the  radio  waves  with  the  ionization 
below  the  reflection  region  causes  the  energy  which  is  given 
to  the  electrons  to  be  lost  in  collisions  with  the  neutral 
atmosphere.  The  effect  of  this  loss  of  energy  is  the 
aUSsorption  of  radio  waves,  and  this  ad>sorption  coefficient, 
which  decreases  with  increasing  frequency,  achieves  such  a 
magnitude  that  propagation  will  obstruct  signal  reception. 

In  the  presence  of  the  earth' s  magnetic  field  and  with  its 
dispersive  properties,  the  ionosphere  becomes  both  anisotropic 
and  doubly  refracting.  So,  in  this  environment,  the 
ionospheric  index  of  refraction,  which  determines  the  radio 
wave  properties  at  any  point  in  the  ionosphere,  is  a  pareuneter 
which  is  a  complex  function  of  four  independent  varied^les  that 
can  be  described  as  functions  of  five  parameters ;  (1)  the 
electron  density  of  the  ionosphere;  (2)  the  atmospheric 
collision  frequency;  (3)  the  radio  wave  angular  frequency;  (4) 
the  ionospheric  angular  gyrofrequency ,  and  (5)  the  angle 
between  the  earth's  magnetic  field  and  the  direction  of  the 
radio  propagation. 

By  using  Snell's  Law,  and  given  the  index  of  refraction  as 
a  function  of  the  radio  wave  position  in  the  ionosphere,  the 
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propagation  of  the  radio  wave  can  be  traced  by  a  form  of  the 
optical  equation.  By  applying  this  law  to  ionospheric 
propagation,  it  can  be  stated  that  when  the  radio  wave 
transverses  an  imaginary  surface  between  two  regions  of  the 
ionosphere  having  different  refractive  indices,  n^  and  nj,  the 
refracted  wave  leaves  the  surface  at  an  anglu  ^2/  which  is 
related  to  the  radio  wave  angle  of  incidence  with  the  surface, 
and  the  indices  of  refraction  as  follows: 

sin ^ 
sin  (4)2)  nl 


The  properties  of  the  ionosphere  are  complex  and  quite 
unstable.  These  properties  affect  the  ray  path  of  the  HF  radio 
wave  in  the  ionosphere,  so  a  complete  description  of  the 
ionosphere  must  take  into  account  this  variability,  part  of 
which  is  cyclic  and  follows  predictaJble  trends,  and  part  of 
which  is  transient.  These  cyclic  variations  are  due  to  the 
earth' s  rotation  (day) ,  to  the  rotation  of  solar  active 
regions  (30  days) ,  to  seasons  (12  months) ,  and  cyclic 
variations  in  solar  activity  (11  years) .  The  noncyclic 

variations  are  due  to  random  changes  in  atmospheric 

properties,  ionization  transport  by  electrical  current  systems 
in  the  ionosphere,  by  large  scale  winds,  and  transient 
radiations  from  the  sun.  The  noncyclic  variations  are 
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unpredictable.  The  only  exceptions  are  some  ionospheric 
effects  produced  by  transient  solar  radiations. 


Among  the  various  anomalous  conditions 
that  have  been  noted,  one  of  the  more 
equatorial  anomaly  [Ref.  14] . 


of  the  ionosphere 
important  is  the 


45 


APPENDIX  C  -  IONOSPHERIC  SOUNDING 


Ionospheric  sounding  is  a  technique  used  for  sounding  the 
ionospheric  propagation  medivun  characteristics.  The  following 
techniques  have  been  developed  for  ionospheric  sounding 
[Ref.  15] : 

Ionospheric  Pulse  Soxinding:  This  concerns  the  measurement 
of  the  linear  unit  impulse  response  function  for  each 
communications  channel.  The  output  of  the  sounder  is  usually 
displayed  in  the  form  of  an  ionograun.  Three  basic  schemes  are 
used: 

1.  Vertical  Incidence  Pulse  Sounding:  The  sounding  pulse 
is  emitter’  vertically  and  the  reflected  returns  from  the 
ionosph  ^re  are  analyzed  at  a  co-located  or  nearby  receiver 
(Fig.  6)  . 

2.  Oblique  Incidence  Bistatic  Pulse  Sounding:  The  sounding 
pulse  is  emitted  either  over  the  actual  communications  path, 
or  over  an  adjacent  path.  This  method  requires  the  transmitter 
and  receiver  to  be  remotely  synchronized  (Figure  7) . 

3.  Oblique  Incidence  Backscatter  Pulse  Sounding:  As  with 
vertical  sounding,  the  transmitter  and  the  receiver  are  co¬ 
located  or  relatively  close,  but  the  received  signals  are 
reflected  obliquely  from  the  ionosphere  and  are  scattered  by 
ground  irrecpularities  (Figure  0)  . 
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Linear  Sweep  Sounding  or  ' Chirpsounding' :  This  consists  of 
sending  a  low  power,  2-30  MHz,  linear  FM/CW  test  signal  over 
the  conununications  path.  The  sounding  signal  is  tracked  by  the 
time-synchronized  receiver.  It  sounds  the  communications  path 
oblicjuely  and  vertically  (Figure  9)  . 

Channel  Sounding:  This  method  uses  probing  transmissions 
on  only  a  limited  number  of  allocated  frequencies  over  the 
radio  link  of  interest.  It  is  more  correctly  described  as 


channel  evaluation. 


Figure  6  -  Vertical  Pulse  Sounding 


Figure  7 


-  Oblique  Bistatic  Pulse  Sounding 
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Figure  8  -  Oblique  Backscatter  Pulse  Sounding 


Figure  9  -  Block  Diagreun  of  Chirpsounding  System 
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APPENDIX  D  -  GLOSSARY 


•  Ionosphere  .  The  ionosphere  extends 

from  50  Km  to  400  km  in  altitude.  At  these  distances  from 
the  surface,  the  air  is  very  thin  and  the  pressure  is  very 
low:  from  2  mm  to  10"®  mm  of  mercury  is  a  typical  range. 
By  contrast,  sea  level  pressure  is  760  mm.  At  such  low 
pressures  as  those  at  the  upper  levels,  the  atoms  of 
atmospheric  gases  are  readily  ionized  by  the  actinic 
ultraviolet  (UV)  light  from  sunlight.  These  ionized  gases 
expel  electrons  which  can  reflect  radio  waves  at  certain 
frequencies  corresponding  to  their  mean  free  paths  of 
travel  before  collision  with  un-ionized  atoms  or  molecules 
[Ref.  14]. 

•  Skywaves  and  groundwaves  .  Radiated  radio  waves 

are  propagated  in  two  main  modes,  by  groundwaves  which  hug 
the  surface  of  the  ground  in  the  troposphere  and  follow 
the  earth' s  curvature  in  this  manner,  and  by  skywaves 
which  are  propagated  in  straight  lines  until  intercepted 
by  the  ionosphere,  which  reacts  in  a  quasi-optical  manner 
(like  a  prism)  by  bending  their  paths  continuously  toward 
the  earth,  up  to  a  limiting  frequency.  This  bending 
results  in  reflection  back  to  the  ground  [Ref.  14] . 

•  Ionospheric  layers  .  The  ionized  region 

consists  of  mixture  of  free  electrons  and  both  positive 
and  negative  ions  distributed  in  a  plasma-like  medium  at 
very  low  pressure.  Only  the  free  electrons  can  move  with 
sufficient  velocity  to  interact  with  high-frequency  radio 
waves .  Although  the  ionization  extends  continuously 
through  the  ionosphere,  there  are  heights  at  which  the 
free  electron  density  reaches  a  maximum  density.  These 
regions  of  maximum  density  are  called  layers  (D,E,F1,F2) 
[Ref.  14]. 

•  Ionization  density  .  The  density  of 

ionization  or  concentration  of  free  electrons  is  quite  low 
in  the  D- layer.  It  is  the  lowest  in  height  and  is 
protected  from  UV  radiation  by  the  higher  E- layer  and  F- 
layer.  However,  it  tends  to  be  ionized  by  solar  X-rays 
which  penetrate  further  than  the  UV  layer  and  make  this 
layer  a  mixture  of  electrons  and  un-ionized  molecules.  As 
the  latter  are  relatively  dense,  they  absorb  most  of  the 
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electrons  excited  by  radio  waves,  and  thus  remove  energy 
from  the  signal.  This  attenuation  occurs  mainly  during  the 
sunlight  hours.  The  E-layer  has  adequate  electron  mean 
free  path  length  to  resonate  with  HF  waves,  and  therefore 
acts  as  a  reflecting  layer.  However,  it  is  very  thin  in 
height  and  at  night  may  lose  most  of  its  ionization  by 
recombination  of  ions  and  electrons.  The  Fl-layer  is  part 
of  the  total  F-layer  complex.  It,  like  the  E-layer,  is 
predominantly  a  daytime  phenomenon,  the  Fl-layer  is  lower 
than  the  F2-layerf  but  merges  with  it  to  some  degree  at 
night.  The  F2-layer  is  the  highest  known  layer  that 
reflects  radio  waves .  It  is  highly  ionized  and  remains  at 
reduced  density,  even  at  nighttime.  Its  height  is  quite 
varisible,  reaching  the  greatest  altitudes  around  local 
noon  [Ref.  14] . 


•  Es . Sporadic  E  is  the  thin 

layer  of  enhanced  ionization  confined  to  the  E  region 
(100-110  Km)  and  is  thought  to  be  a  result  of  an  enhanced 
concentration  of  ions  caused  by  wind  shears  in  that  region 
[Ref.  6] . 

•  Spread  F  .  Spread  F  is  caused  by 

the  scattering  of  the  signal  from  irregularities  embedded 
in  the  i .nosphere  both  in  depth  and  away  from  the  zenith. 
The  echo  pulse  reflected  from  the  F2  layer  has  a  much 
longer  duration  than  the  transmitted  pulse.  The  echoes 
obtained  from  that  layer  by  vertical  sounding  become 
diffuse  and  of  indefinite  height.  The  effect  in  the  F2 
layer  soon  after  local  sunset,  and  which  appears  to  last 
for  a  few  hours  thereafter  in  a  region  near  the  earth' s 
magnetic  equator  is  known  as  "equatorial  spread  F" 
[Ref.  16]  . 

•  Tilts . Ionospheric  tilts  refer 

to  generally  large-scale  horizontal  gradients  in  the 
electron  density,  such  that  contours  of  constant  electron 
density  are  no  longer  parallel  to  the  earth' s  surface 
[Ref.  6] . 

•  TIDs  .  Traveling  ionospheric 

disturbances  are  essentially  an  ionospheric  manifestation 
of  an  entire  spectrum  (not  necessarily  continuous)  of 
waves  propagating  through  the  atmosphere  [Ref.  6] . 

•  Vertical  sounding  .  Transmission  of  a  HF 

wave  vertically  to  the  ionosphere  is  a  means  of 
determining  its  height  and  frequency-reflecting 
characteristics.  The  treuisit  time  delay  is  a  function  of 
the  height  to  the  maximum  intensity  layer.  The  upper 
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frequency  limit  is  a  measure  of  the  degree  of  ionization 
[Ref.  14]. 

•  Oblique  sounding  .  In  oblique  sounding, 

the  angle  of  incidence  of  the  transmitted  radio  wave  is 
less  than  90  degrees.  Thus,  the  point  of  return  to  ground 
is  remote  from  the  transmitter.  This  wave  can  then  be 
reflected  alternately  between  the  earth  and  the  ionosphere 
several  times.  Through  variations  in  the  ionosphere  and/or 
multi-angle  radiation  from  the  antenna,  many  reflection 
paths  can  be  taken  by  the  reflected  waves.  It  is  also 
possible  for  individual  reflections  from  the  E-layer  and 
F-layer  to  interact.  These  can,  and  often  do,  exist 
simultaneously,  causing  multisignal  reflection  from  a 
single  transmitter  at  a  remote  receiving  point.  The 
oblique  type  of  sounding  is  valuable  for  the  investigation 
of  propagation  paths  over  any  practical  path  using  the 
ionosphere  [Ref.  14]  . 

•  Virtual  height . Normally  designated  h', 

the  virtual  height  is  the  height  at  which  a  signal 
transmitted  toward  the  ionosphere  appears  to  be  reflected. 
It  is  also  the  height  obtained  by  vertical  incidence 
sounders,  since  the  time  delay  of  the  actual  path  of  the 
signal  is  nearly  the  same  as  the  apparent  path  of  the 
signal,  provided  that  the  ionosphere  is  treated  as  a  plane 
mirror  at  the  apex  of  the  apparent  path.  As  the  frequency 
of  a  vertically  incident  transmitted  signal  is  increased, 
the  measured  virtual  height  of  the  signal's  reflection 
increases.  Although  the  actual  speed  in  the  ionosphere  is 
less  than  in  free  space;  the  height  is  calculated  as  if 
the  medium  were  free  space,  hence  the  adjective  "virtual." 
[Ref.  14] 

•  hmF2,  hmFl,  and  hmE  .  The  height  of  the 

maximum  electron  density  in  the  F2  and  FI  layers  and  the 
E-region  respectively. 

•  Critical  frequency  .  It  is  the  highest 

frequency  which  is  reflected  back  to  ground  during 
vertical  sounding.  At  oblique  angles  of  incidence,  with 
longer  paths  traversed  through  the  ionosphere,  higher 
frequencies  can  be  reflected  back  to  earth,  and  calculated 
as  a  product  of  the  critical  frequency  with  the  sine  of 
the  oblique  angle  of  incidence  [Ref.  14] . 

•  Maximum  us2d)le  frequency  (MUF)  . . .  The  maximum  frequency 
that  can  be  used  for  long  range  communications  (>3000  km) 
can  be  higher  than  the  critical  frecjuency  of  vertical 
sounding.  This  maximum  useible  frequency  for  oblique 
incidence  transmission  is  calculated  as  the  product  of  the 
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critical  frequency  and  the  secant  of  the  oblique  angle  of 
incidence  (secant  law)  .  This  formula  assvunes  a  flat  earth 
and  a  similarly  flat  F-layer.  The  curvature  of  the  earth 
introduces  serious  errors  beyond  a  few  kilometers. 
Frequencies  between  the  highest  and  lowest  usedjle  limits 
will  be  transmitted  and  received.  However,  many 
frequencies  in  this  band  are  subject  to  multiple 
reflections  which  degrade  the  quality  of  reception.  The 
shortest  distance  at  which  any  frequency  can  be  received 
by  reflections  is  the  "skip"  distance.  The  skywave  signals 
cannot  be  received  at  shorter  distances  than  the  "skip" 
distance  [Ref.  14] . 

•  Multipath  interference  .  When  the  frequency  of 

transmission  is  appreciably  less  than  the  MUF,  the  signals 
are  received  by  several  propagation  paths  having  different 
delay  times.  These  time  differences  result  in  distortion 
of  phase  of  the  signal  modulation,  which  can  be  severe  in 
the  case  of  navigational  signals,  and  cause  selective 
fading  of  side  frequencies  in  the  modulation  sidebands. 
Since  the  propagation  paths  are  very  sensitive  to 
geometrical  variations  of  the  reflection  angles, 
turbulence  in  the  ionosphere,  which  is  common,  causes 
either  rapid  or  slow  variations  of  flutter-fading  through 
changing  delay  times .  While  frequencies  near  the  MUF  are 
less  susceptible  to  multipath  transmission,  they  are 
sxabject  to  flutter- fading  due  to  temporary  variations  in 
the  locally  ionized  regions  [Ref.  14]  . 

•  Ordinary  and  Extraordinary  Rays  .  .  The  permanent  magnetic 

field  of  the  Earth  makes  the  ionized  gas  of  the  ionosphere 
an  anisotropic  medium  (different  properties  in  different 
directions) .  When  a  linearly  polarized  wave  is  propagated 
in  a  arbitrary  direction  relative  to  the  direction  of  the 
magnetic  field  of  the  earth,  the  wave  is  split  into  two 
elliptic  rays,  the  ordinary  (O)  and  the  extraordinary  (E) 
rays.  The  longitudinal  electric  field  components  of  both 
the  E  and  0  rays  are  in  phase  with  the  transverse 
component,  oriented  in  the  direction  of  the  magnetic  field 
of  the  earth.  These  rays  are  propagated  with  different 
velocities  and  paths,  and  the  major  axes  of  the  ellipses 
of  polarization  are  at  the  right  angles  to  each  other, 
with  the  resultant  vectors  rotating  in  opposite  directions 
[Ref.  17]  . 

•  Gyro- frequency  .  The  rotation  of  the 

electrons  under  the  effect  of  the  terrestrial  magnetic 
field  exhibits  resonance  characteristics.  This  is  why  it 
is  sometimes  referred  to  as  gyromagnetic  resonance,  and 
the  respective  frequency  is  called  the  gyro- frequency .  The 
average  value  of  the  constant  magnetic  field  of  the  earth. 
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is  40  A/m  and  gives  1.4  MHz  as  the  gyro- frequency .  It 
may  be  expected  that  waves  at  a  frequency  of  eJaout  1 . 4  MHz 
will  experience  increased  attenuation  in  the  ionosphere 
due  to  collision  losses  [Ref.  17] . 

•  foF2,  foFl,  and  foE  .  The  critical  frequency 

of  the  ordinary  wave  for  the  F2  and  FI  layers  and  the  E 
region,  respectively  [Ref.  18] . 

•  M(n)X  .  The  maximum  usable 

factor  for  a  path  of  "n"  kilometers  for  transmission  by 
the  "X"  layer.  For  exaunple,  M(3000)F2  represents  the 
maximum  usable  factor  for  a  path  of  3000  km  for 
transmission  by  the  F2  layer  [Ref.  18] . 

•  Geomagnetic  Latitude  .  The  boundaries  between 

the  low,  middle,  and  high  latitudes  are  not  fixed,  but 
vary  with  local  time,  season  of  year,  and  solar  cycle.  So, 
the  "low  latitude"  can  be  viewed  as  extending  from  between 
+/-  30°  geomagnetic  latitude,  the  "high  latitude"  greater 
than  60°  geomagnetic  latitude,  and  "middle  latitude"  as 
the  geomagnetic  latitudes  between  30°  and  60°.  Because  the 
geomagnetic  and  geographic  poles  do  not  coincide,  there 
are  substantial  differences  between  the  geographic  and 
geomagnetic  latitudes,  particularly  in  the  Americas  and 
throughout  Asia. 

•  Covariance . Covariance  measures  how 

strongly  two  random  variables  X  and  Y  are  related  to  one 
another  [Ref.  4] . 
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APPENDIX  E  -  IONOSPHERIC  DATA  RESEARCH 


In  the  region  near  the  earth's  magnetic  equator,  several 
magnetic  and  ionospheric  phenomena  occur  that  appear  to  be 
peculiar  to  that  region.  Evidently,  much  investigation  must  be 
done  before  these  occurrences  become  fully  explained. 

"The  Condor  Equatorial  Spread  F  Campaign, "  published  at 
the  Journal  of  Geophysical  Research,  Vol.  91,  No  A5,  (May  1, 
1986) ,  summarizes  the  main  results  of  many  experiments 
designed  to  study  the  South  American  equatorial  F  region 
[Ref.  19]. 

Based  on  that  piablication  the  following  agencies,  were 
contacted  in  June,  1990  to  suggest  a  source  of  vertical 
incidence  ionosonde  data  or  incoherent  scatter  measurements 
which  could  be  use  to  ested^lish  the  normal  vs.  tilted 
ionospheric  profiles  in  that  region.  Unfortunately,  none  of 
these  requests  produced  positive  results. 

•  Air  Force  Cambridge  Research  Ledsoratories,  Bedford, 
Massachusetts 

•  Center  for  Atmospheric  and  Space  Sciences,  Utah  State 
University,  Logan,  Utah 

•  Emmanuel  College,  University  of  Massachusetts,  Boston, 
Massachusetts 
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•  Institute  Geofisico  del  Peru,  Lima,  Peru^ 

•  Ionospheric  Physics  Division,  Air  Force  Geophysics 
Laboratory,  Hanscon  Air  Force  Base,  Massachusetts 

•  Los  Alaunos  National  Laboratory,  Los  Alcunos,  New  Mexico 

•  National  Space  Research  Institute  (INPE) ,  Sao  Jose  dos 
Ccunpos,  SP,  Brazil  ^ 

•  School  of  Electrical  Engineering,  Cornell  University, 
Ithaca,  New  York 

•  Space  Physics  Research  Laboratory,  University  of  Michigan, 
Ann  Arbor,  Michigan 

•  The  University  of  Texas  at  Dallas,  Dallas,  Texas 


^  Dr.  Ronald  F.  Woodman,  Director  of  the  "Institute  Geofisico 
Del  Peru"  reported  that  data  from  "Radio  Observatorio  de  Jiceunarca" 
are  single  station  (no  tilt  information) . 

^INPE  was  recommended  by  Mr  Jorgen  Buchau  from  the  Ionospheric 
Physics  Division,  Hanscon  Air  Force  Base,  Massachusetts. 
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